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SUPPLEMENTAL MATERIAL 
THz absorption spectra in the two Ba3NbFeaSi20i4 crystals 

Figure 1 presents the absorption/absorbance spectra (left/right scales) measured with the 6 BMS on two 
different samples. The absolute absorbances, ad, with a the absorption and d the travelled distance through the 
sample, were determined by measuring the transmission through a 2 mm diaphragm as a reference and the sample 
transmission through that same diaphragm. Interference pattern from the sample surfaces is present for sample 2 
(lower panel). For further analysis, this interference pattern was removed using a band block FT filter. The larger 
oscillations observed below 30 cm -1 for sample 1 (upper panel) are instrumental. For further analysis, they were 
removed using the sample at a different temperature as a reference. The stronger absorption observed for e//c is 
attributed to a (multi-)phonon band. 
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FIG. 1: Absorption spectra in Ba3NbFe3Si20i4 measured with the 6 /im BMS. The temperature was varied from 6 
K to 300 K. Three different polarizations of the THz electromagnetic field as regards the crystal c-axis were measured. Upper 
panel: measurements on sample 1, (e//c h_Lc) and (e±c h//c). Lower panel: measurements on sample 2 (e±c h±c). 
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Interaction between R z and the external fields in Eq.l 



In isotropic or cubic systems, the external field conjugated to rotations is a torque. However, in systems with 
reduced symmetry, the linear coupling between rotations R and linear displacements can be allowed. This is precisely 
our case, since our system initially has D% point-group symmetry which is further reduced to G%. This means that R 
in our system is not a pure rotation anymore, so that the standard forces acting on the above displacements also act 
on R. Polar displacements implying a relative motion of charges are among the displacements linearly coupled to R. 
Thus, the electromagnetic radiation is coupled to rotations via this motion of charges, which gives rise to the unusual 
electromagnetic activity of the atomic rotations in our system. 

More specifically, the charges feel the electromagnetic radiation through the Lorentz force: f = — | ^ — eS7<f> + 
-fx (V x A), where r represents the position of the charge e, while <j> and A are the scalar and vector potentials 
of the radiation respectively [1]. In the Weyl gauge (<fi = 0), the vector potential of a monochromatic plane wave is 
A(r, t) = a e^ k r ~"*) ~ (1 + ik • r)Ao(t) in the long-wavelength limit. Here k = nto/c and A (t) = aoe -1 "*, with n 
being the refractive index of the medium and ao a constant complex vector. In this limit, the first term in the vector 
potential generates the electric field of the wave E = i(k/n)Ao, while the magnetic field H = ik x Ao is obtained 
from the second one. Keeping this in mind, the vector potential can be rewritten as A ~ + H x r + n£(E • r). 
Accordingly, the Lorentz force becomes f ~ eE + |r x H + zek(E • r) + -f x H. We note that the second and the last 
contribution to this force are sort of dual contributions. Interestingly, the last term was invoked by Dzyaloshinskii 
and Mills in order to explain the unexpected paramagnetic behavior observed in nonmagnetic ferroclcctrics [2]. In 
our case, the second term is the key one for the physical interpretation of the absorption spectrum of our system. In 
fact, in our measurements, the last two contributions to the Lorentz force produce subdominant effects and therefore 
can be ignored. Furthermore, this force acts directly on the charge distribution of the system, and therefore on its 
electric polarization P = J~l e^. But since this polarization is linearly coupled to atomic rotations in our system, this 
gives rise to an effective force acting on the rotations as explained above. The force acting the rotations R z about the 
z-axis, in particular, can be written as 



Here a represents the strength of the coupling between R z and P z , up z is the frequency (or stiffness) associated to 
P z , and g is a gyromagnetic factor. This force can be obtained as / c ft = ^ nt , where Lj nt is the Lagrangian in Eq. 
(1) of the main text. 

If is uniform, both electric and magnetic-field components of the THz radiation will couple to q = rotations 
due to momentum conservation. However, if P± is modulated with a finite wavevector Q, then L lnt becomes 



in Fourier space, where a\ = a/u) P (0) and ot2 = ga/cup (Q). In this case, the electric field couples to q = rotations 
while the magnetic field does to q = ±Q ones, thus probing different resonant energies. 
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— E z +g— (P± x H_l) z . 




(1) 



lira ai E z (k)R z (-k) + a 2 (P±(-Q) x H±(k)) x R,(Q - k) + (Q -> -Q) 



(2) 
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We have determined the terahertz spectrum of the chiral langasite Ba3NbFe3Si20i4 by means of 
synchrotron-radiation measurements. Two excitations are revealed that are shown to have a different 
nature. The first one, purely magnetic, is observed at low temperature in the magnetically ordered 
phase and is assigned to a magnon. The second one persits far into the paramagnetic phase and 
exhibits both an electric and a magnetic activity at slightly different energies. This magnetoelectric 
excitation is interpreted in terms of atomic rotations and requires a helical electric polarization. 
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The electric-field control of spins and the converse 
magnetic-field control of electric dipoles inspire a num- 
ber of hybrid technologies and motivates fundamental 
research on multiferroics and magnetoelectric materi- 
als [1]. These magnetoelectric couplings produce strik- 
ing phenomena both at the static and at the dynam- 
ical level. A prominent example is the electric-charge 
dressing of magnons, resulting in the so-called clcctro- 
magnons, that has been demonstrated experimentally in 
multiferroics [2, 3]. This dressing enables the electric- 
field control of spin-waves and produces enhanced op- 
tical responses with potential applications in photonics 
and magnonics [4]. Here we report experimental evi- 
dence of the dual phenomenon, that are atomic vibra- 
tions dressed with currents, hence magnetically active, 
in the terahertz (THz) spectrum of the chiral compound 
Ba3NbFe3Si20i4. Based on symmetry arguments, wc 
propose that this new type of magnetoelectric excita- 
tions is associated to structural rotations; their unex- 
pected magnetoelectric activity is triggered by a sponta- 
neous helical polarization in the paramagnetic phase of 
the system. 

When chirality meets magnetic order, unusual mag- 
netic and electric properties often appear encouraging 
the search for novel multifcrroic and magnetoelectric phe- 
nomena. The highest magnetically-induced polarization 
observed so far has been reported in CaMn70i2 [5] where, 
interestingly, axial distortions of the lattice play a key 
role in generating ferroclectricity out of the chiral mag- 
netic structure of the system [6, 7]. We hereafter focus on 
the Ba3NbFe3Si20i4 langasite, which crystallizes within 
the non-centrosymmetric space group P321 [8] thus dis- 
playing structural chirality (see Fig. 1). The remark- 
able magnetic properties of this system are due to the 
Fe 3+ ions. These form triangles stacked along the c-axis 



and compose a triangular lattice in the (a, b) plane (see 
Fig. 1). The magnetic exchanges along the c-axis are 
mediated by twisted Fe-O-O-Fe paths. Below the Neel 
temperature Tn = 27 K, the system develops a chiral 
magnetic order with a 120° spin structure within the Fe- 
triangles that rotates along the c-axis and forms magnetic 
helices [9]. In addition, below T^r, the system shows a 
static magnetoelectric effect [10] and chiral spin- wave ex- 
citations [11-13]. A weak electric polarization along the 
c-axis has also been reported [14], although our own mea- 
surements show no uniform polarization cither along the 
c-axis nor perpendicular to it. 

We have studied Ba 3 NbFc3Si20i4 by means of 
synchrotron-radiation measurements in the THz regime. 
This technique is particularly well suited to investigate 
magnetoelectric excitations [15] since (i) the THz energy 
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FIG. 1: (a) Ba3NbFe3Si2 0i4 langasite crystallographic struc- 
ture, (b) Relative orientation of the electric e and magnetic 
h fields of the THz wave with respect to the crystal c-axis for 
samples 1 and 2. 
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range coincides with that of magnons and phonons and 
(ii) the THz-wave delivered by the synchrotron radiation 
(98 % linearly polarized) carries both an electric and a 
magnetic field, so that rotating the crystal with respect 
to the wave polarization direction allows the unambigu- 
ous determination of the excitation fields. We have used 
these selection rules and the temperature dependence of 
the excitations to elucidate their nature and the mech- 
anisms at play in the dynamics. This is further ascer- 
tained by comparison with inelastic neutron scattering 
measurements. 

Two plaqucttes of mm size were cut from the same sin- 
gle crystal [9] with a thickness of 450 fim and 300 /xm; the 
c-axis was oriented either in plane or out of plane in order 
to probe all possible geometries of the THz electric and 
magnetic fields as regards the crystal c-axis: (e//c h_Lc) 
and (e_Lc h//c) for one crystal, (e_Lc h_Lc) for the 
other crystal (see Fig. 1). THz absorption spectra were 
obtained in the energy range 8-60 cm" 1 (0.24-1.8 THz) 
by measuring the transmission at the AILES beamline 
of Synchrotron SOLEIL [16]. A Bruker IFS125 interfer- 
ometer equipped with a pulse tube cryostat was used, 
combined with a Helium pumped bolometer. The 10- 
60 cm -1 energy range was explored at a resolution of 
0.5 cm -1 using a 6 fim thick silicon-mylar multilayered 
beamsplitter (BMS). The 8-18 cm -1 energy range was 
explored with more sensitivity by using a 125 fim thick 
mylar BMS. 

The THz spectra measured at various temperatures 
from 300 K to 6 K (for more details, see Supplemen- 
tal Material) reveal an absorption roughly three times 
larger for e//c than for e±c, meaning that the dielec- 
tric response of the material is considerably stronger 
for this direction of the THz electric field. This addi- 
tional absorption arises from (multi-)phonons excited by 
e//c, that contribute as a broad background centered 
at 40 cm" 1 growing with increasing temperature, as ex- 
pected for thermally populated phonons. 

Since the magnetic contribution is weak compared to 
the instrument and phononic background, two different 
methods were used to enhance the magnetoelectric re- 
sponse of the samples. For the 6 fim BMS measure- 
ments, we used the same procedure as in [17]: we es- 
timated the extrapolated absorption at a given tempera- 
ture To, a(To), as the common signal to the a(To)-a(T) 
curves for all measurements temperatures T. A linear 
background was then subtracted. On the other hand, for 
the 125 fim BMS measurements, the sample at 45 K was 
used as a common reference. Interference patterns aris- 
ing from coherent reflections on the disk surfaces, when 
present, were removed using a band block FFT filter. 
The extrapolated absorption for all the different combi- 
nations of the THz electric field e and magnetic field h 
alignment is displayed on Fig. 2. 

The most stiking feature is visible for (e//c h_Lc): 
a large asymmetric absorption is clearly observed that 
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FIG. 2: THz excitations probed with the 6 fim BMS. 

(a) Extrapolated absorption spectra for (e//c h_Lc) at dif- 
ferent temperatures. Insert : gaussian fits. (b)(c)(d) Tem- 
perature dependance of the fits position, width and area. 
(e)(f) Same as (a) for the other THz polarizations. 
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FIG. 3: THz excitations probed with the 125 fi 
BMS. (a) Relative absorption spectra, abs(T)-abs(45 K) 
(e_Lc h//c) at different temperatures. Insert : gaussian fits. 
(b)(c)(d) Temperature dependance of the fits position, width 
and area. 



results from two gaussian contributions (see insert of 
Fig. 2a). At 16 K for instance, the two peaks are found at 
23 and 29 cm" 1 with a respective width of 10 and 5 cm -1 . 
Rotating the electric field perpendicular to c leaves with 
a single broad absorption centered around 25 cm" 1 (see 
Fig. 2e). When an additional rotation is performed and 
the magnetic field lies parallel to c no clear absorption 
is observed (see Fig. 2f). These measurements evidence 
two excitations quite close in energy, El and E2, that 
persist in the paramagnetic phase up to 100 K, that is 
almost four times higher than T^r (see Fig. 2). This is a 
strong indication that they have a common origin which 
is not related to the long range magnetic order although 
they lie in the magnon energy range. It leads us to in- 
troduce a new characteristic temperature in the system, 
T p ~ 100 K. 

The overall behavior of these excitations presents some 
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discrepancies: the energy of excitation El remains con- 
stant while for E2 it increases at high temperature. The 
life time (inverse of the absorption peak width) starts 
to decrease above Tjv for excitation El while it remains 
constant for E2. The most remarkable difference con- 
cerns their selection rules: El is activated by the mag- 
netic field h_Lc, E2 by the electric field e||c [see Fig 2 
panels (a), (e) and (f)]. Note that both El and E2 are 
clearly seen in panel (a), El is still visible in panel (e), 
and no feature is detectable in panel (f). 

Another weaker and even sharper excitation (E3) is 
found at lower energy for (e_Lc h//c), using the more 
sensitive 125 /im BMS, resulting from two close gaussian 
contributions. At 10 K, they are observed at 12.9 cm -1 
and 14.1 cm -1 , both with a width of 0.8 cm -1 (sec 
Fig. 3a). This excitation shifts to lower energy and 
broadens as the temperature is increased to finally disap- 
pears at 24 K, just below Tn- Contrary to the previous 
excitations, the temperature dependence of this excita- 
tion is characteristics of a magnon. 

This was confirmed from the comparison with the 
magnon spectra obtained through inelastic neutron scat- 
tering [11-13, 18]. The THz measurements probe the 
first Brillouin zone center (0,0,0) where no neutron mea- 
surements are available. We therefore compared them 
to neutron spectra in other zone centers. Selected mea- 
surements performed on the time-of-flight IN5 and on 
the triple-axis IN12 spectrometers at the Institut Lauc 
Langevin are shown in Fig. 4. Two main branches emerge 
from the magnetic satellites at +t with r=l/7 [19]. The 
lower energy branch is associated with the dynamical cor- 
relation function of spin components along the c-axis. 
The upper energy branch is associated with those in 
the triangle plane. An intense signal is observed around 
13 cm -1 at the zone center, mainly from the lower branch 
(see Fig. 4b). The tail of this signal and the less in- 
tense neutron signal up to 25 cm -1 are attributed to the 
upper branch. The first neutron peak at 13 cm -1 per- 
fectly agrees with the weak excitation E3 observed in the 
THz measurements in terms of energy and spin compo- 
nent involved (along the c-axis, hence excitable by h//c). 
The other magnonic signals (involving spin components 
in the (a,b) plane) are not clearly distinguished in the 
THz measurements for h±c, presumably because the 10 
times stronger excitations El and E2 mask everything 
(see Fig. 2a and Fig. 2e). These remarkable larger THz- 
excitations El and E2 are clearly absent from the neutron 
data of Fig. 4 and also from references [11, 12, 18]. 

In multiferroics, magnons can be electro-active at par- 
ticular points of the Brillouin zone (related to the spiral 
wave vector or at the Brillouin zone edges) [2, 3, 20] and 
these so called electromagnons can persist in the non- 
multiferroic, yet magnetically ordered phases of these 
compounds [2, 3, 21]. In contrast, for the langasite com- 
pound, the THz excitations El and E2 persist up to 
100 K 4T/v when no well-defined magnetic excitation 
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FIG. 4: Magnons probed by inelastic neutron scatter- 
ing, (a) Magnons dispersion along the (0,-1, -£) direction in 
the reciprocal space. The black line indicates the Brillouin 
zone center, (b) Energy scan at the reciprocal lattice zone 
center (0,1,0). The colored dash lines give the energy posi- 
tion of excitations El, E2 and E3 seen in THz measurements 
while the colored rectangles indicate their width. 

exists anymore [22]. Thus, whereas the THz-excitation 
(E3) at 13 cm -1 can be unambiguously identified as a 
conventional magnon, the magnonic origin of the THz- 
excitations El and E2 can be discarded. 

These excitations are more likely connected to the lat- 
tice [23], although their relatively low energies and the 
unusual behavior described above is different from that 
of standard phonons. The natural candidates to explain 
these features are axial distortions of the lattice, whose 
role has been recently highlighted in this class of mate- 
rials [6, 24, 25]. These are rotations R z about the c-axis 
whose interaction with the electric and magnetic fields 
of the external radiation is described by the effective La- 
grangian: 

L int = -^-[e z +g(P ± xh ± ) z }R z . (1) 

This interaction results from the linear coupling between 
rotations and polar displacements allowed by the P321 
space-group symmetry of our system (see Supplemental 
Material for details on the derivation of this Lagrangian) . 
Here a represents the strength of this coupling, P is the 
electric polarization, with wp z being its characteristic fre- 
quency (or energy), and g is a gyromagnctic factor. 

The excitation of R z by the electric-field component of 
the radiation observed in our THz measurements is de- 
scribed by the first term in Eq. (1). This follows imme- 
diately from the aforementioned linear coupling between 
R z and P z . The additional response of R z to the mag- 
netic field is due to the second term in Eq. (1). This re- 
sponse, in contrast, is far more subtle. It requires a finite, 
static polarization perpendicular to the c-axis, P± ^= 0, 
to be operative. At room temperature, no perpendicu- 
lar polarization is allowed by the symmetry. However, at 
low temperatures, recent neutron scattering results are 
not compatible with a perfect 120° spin structure, im- 
plying the loss of the three- fold symmetry axis [22]. This 
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P321 — ► C2 symmetry breaking implies the appearance 
of a spontaneous polarization perpendicular to the c-axis 
that, in turn, activates the coupling of R z with the 
magnetic field in Eq. (1). This is in perfect agreement 
with the gradual appearance of the magneto-active ex- 
citation El that we observe at 24 cm -1 below 100 K 
(see Fig. 2) . The concomitant appearance of the electro- 
active excitation E2 at 100 K can be explained as due 
to the increase of the overall effective coupling a/tu P at 
the onset of the new phase. This increase is likely a by- 
product of the softening of uip z due to its dependence on 
P±. In fact, we observe a similar softening in the energy 
of the electro-active excitation itself that is in tune with 
this hypothesis (see Fig. 2b). 

Finally, we note that in our spectra the energy of the 
electro-active and magneto-active excitations are slightly 
different. This difference suggests that the spontaneous 
polarization that appears below 100 K is modulated with 
a finite wavevector Q. Thus, according to Eq. (1), the 
electric field couples to q = uniform rotations while the 
magnetic field does to q = ±Q rotations whose resonant 
energy is lower (see Fig. 5). The detailed Ginzburg- 
Landau analysis of the P321 — > (72 transition confirms 
this possibility. The perpendicular polarization P± can 
be seen as the order parameter driving the P321 — > C2 
transition. The Ginzburg-Landau free energy then can 
be written in terms of (P x ,P y ) = (p cos 8 , p sin 9) as: 

F=^p 2 + 7 p 3 cos36+^p 4 

-\p 2 (d z 8)+ C -[(S7p) 2 +p 2 (V8) 2 ]. (2) 

Here, as customary, a is assumed to be the only coeffi- 
cient that depends on the temperature and changes sign 
as a = a'(T — T ), where T can be seen as the nom- 
inal transition temperature. If the space variations of 
the order parameter are ignored, then the cubic term in 
the free energy implies rather a first-order transition at 
Tq = T a + 2"/ 2 /(a'b). The presence of the Lifshitz invari- 
ant —Xp 2 (d z 0), however, leads to a second-order transi- 
tion at higher temperature, which sets up helical spon- 
taneous polarization (P x ,P y ) = (P cos Qz, Pq sinQz) at 
T P = T Q + X 2 /(a'c) (see Fig. 5). Here P 2 = \T-T P \/(a'b) 
and Q = A/c. This will actually be the case if 2j 2 jb < 
A 2 /c (that is, if the expected discontinuity at the first- 
order transition is relatively weak). This type of helical 
polarization has been discussed theoretically in the con- 
text of incommensurate phases, although no experimen- 
tal realization has been reported so far [26]. Here Tp 
is different from the direct P321 — > C2 transition tem- 
perature, which is avoided due to the formation of the 
intermediate helical state. The C2 symmetry is obtained 
when the helix becomes commensurate. This is expected 
to take place at a lower temperature, which may well 
coincide with the magnetic transition temperature. 
The essential ingredients that give rise to the mag- 




TEMPERATURE 

FIG. 5: Temperature evolution of the electric and 
magnetic phases. From top to bottom: symmetry, order 
parameter, and corresponding excitations. From right to left 
(decreasing temperature): T > T p paramagnetic and para- 
electric state, T p > T > Tjv helicoidal state of polarization, 
Tjv > T additional helicoidal magnetic state. The P321 to 
C2 symmetry breaking is schematized by the loss of the C3 
axis in the Fe 3+ triangles (from equilateral to isosceles). The 
curvature of the dispersion curve E Tot (q) is deduced from the 
relative energies of the magnetic and electric active peaks El 
and E2. 

netoelectric excitations in our system are also present 
in the newly reported ferroaxial multifcrroic CaMn70i2 
and related compounds [5-7, 25]. There, the linear cou- 
pling between the electric polarization and macroscopic 
structural rotations [6, 25], is preceded by a structural he- 
lix due to the incommensurate modulation of the atomic 
bonds [7], in CaMn70i2 at least. However a chiral mag- 
netic order is necessary to activate the polar-axial cou- 
pling via the trilinear invariant aP z R z , where a repre- 
sents the chirality of the magnetic structure. In the lan- 
gasites, the lattice itself is chiral, so that the polar-axial 
coupling is "switched on" already in the paramagnetic 
phase. Thus, while structural rotations in both these 
systems can be expected to have similar magnetoelectric 
responses, they can reveal different key aspects associ- 
ated to structural vs. magnetic chirality. 

In summary, our THz investigation of the langasite dy- 
namical properties has revealed a new kind of excitation, 
associated to atomic vibrations, which is both electric 
and magnetic active, at slighly different energies. These 
findings demonstrate that not only magnons but also 
atomic vibrations can acquire a magnetoelectric charac- 
ter, thus opening new routes to carry and process infor- 
mation. This striking magnetoelectric activity has been 
related to the bilinear coupling between the atomic rota- 
tions allowed by the chiral symmetry of our system and 
a perpendicular electric polarization. The predicted heli- 
cal polar order may itself exhibit original properties such 
as dielectric chirality and sets off new magnetoelectric 
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functionalities that can be common among symmetry- 
equivalent systems. 

We are very grateful to J. Ollivier, M. Enderle and P. 
Steffens for their help during the neutron experiments on 
IN5 and IN12. We thanks P. Bordet for discussions on 
the langasite structure. We acknowledge J. Debray, J. 
Balay and A. Hadj-Azzem for the crystals preparation. 



* Electronic address: sophie.debrion@grenoble.cnrs.fr 
t Institut Laue Langevin, 6 rue Jules Horowitz, BP 156, 
F-38042 Grenoble Cedex 9, France 
[1] For reviews, see e.g. M. Bibes, and A. Barthelemy, IEEE 
Trans. Electron. Devices 54, 1003 (2007); S.-W. Cheong 
and M. Mostovoy, Nature Materials 6, 13 (2007); W. 
Eerenstein, N. D. Mathur, and J. F. Scott, Nature 442, 
759 (2006); M. Fiebig, J. Phys. D: Appl. Phys. 38, 
R123 (2005); N. A. Spaldin and M. Fiebig, Science 309, 
391(2005). 

[2] A. Pimenov et al., Nature Physics 2, 97 (2006). 

[3] For reviews, see e.g. N. Kida et al., Soc. Am. B 26, A35 

(2009); A. M. Shuvaev, A. A. Mukhin, and A. Pimenov, 

J. Phys.: Condens. Matter 23, 113201 (2011); M. Ken- 

zelmann, Physics 4, 88 (2011). 
[4] P. Rovillain et al., Nature Materials 9, 975 (2010); I. 

Kezsmarki et al., Phys. Rev. Lett. 106, 057403 (2011); 

Y. Takahashi et al., Nature Physics 8, 121 (2012); S. 

Bordacs et al., Nature Physics 8, 734 (2012). 
[5] G. Zhang et al., Phys. Rev. B 84, 174413 (2011). 
[6] R.D. Johnson et al., Phys. Rev. Lett. 108, 067201 (2012). 
[7] N.J. Perks et al., Nature Comm. 3, 1277 (2012). 



[8] B. V. Mill' et al., Russian Journal of Inorganic Chemistry, 

43, 1270 (1998). 
[9] K. Marty et al., Phys. Rev. Lett. 101, 247201 (2008). 
[10] K. Marty et al., Phys. Rev. B 81, 054416 (2010). 
[11] M. Loire et al., Phys. Rev. Lett. 106, 207201 (2011). 
[12] J. Jensen, Phys. Rev. B 84, 104405 (2011). 
[13] V. Simonet, M. Loire, and R. Ballou, Eur. Phys. J. Spe- 
cial Topics 213, 5 (2012). 
[14] H. D. Zhou et al., Chemistry of Materials 21, 156 (2009). 
[15] S.-I. Kimura, H. Okamura.,arXiv: 1210.5310. 
[16] P. Roy et al., Infrared Physics and Technology 49, 139 
(2006). 

[17] S. de Brion et al., Phys. Rev. B 75, 094402 (2007). 

[18] C. Stock et al, Phys. Rev. B 83, 104426 (2011). 

[19] The — t satellites in the presented Brillouin zones are 

much weaker as a consequence of the structural chirality 

[11]. 

[20] H. Katsura, A. V. Balatsky, and N. Nagaosa, Phys. Rev. 

Lett. 98, 027203 (2007); A. Cano and E. I. Kats, Phys. 

Rev. B 78, 012104 (2008); R. Valdes Aguilar et al., Phys. 

Rev. Lett. 102, 047203 (2009). 
[21] A. Cano, Phys. Rev. B 80, 180416 (2009); S. Seki et al., 

Phys. Rev. Lett. 105, 097207 (2010). 
[22] L. Chaix et al, to be published. 

[23] Note that the neutron spectra measured only within 
the first and second Brillouin zones can scarcely probe 
phonons. 

[24] M. Fiebig, Nature Materials 10, 339, (2011); P. Ghosez 
and J.-M. Triscone, Nature Materials 10, 269 (2011). 

[25] A.J. Hearmon et al., Phys. Rev. Lett. 108, 237201 (2012). 

[26] A. P. Levanyuk, in Incommensurate phases in dielectrics, 
vol. 1, Eds. A. P. Levanyuk and R. Blinc (Amsterdam, 
North- Holland, 1986). 



